m 



^ 



X 

5-H 



IPHT 
LPTENS 10/23 



A matrix model for the topological string II 

The spectral curve and mirror geometry 



O B. Eynard 1 , A. Kashani-Poor 2,3 , O. Marchal 1,4 

1 Institut de Physique Theorique, 

CEA, IPhT, F-91191 Gif-sur-Yvette, France, 

CNRS, URA 2306, F-91191 Gif-sur-Yvette, France. 

2 Institut des Hautes Etudes Scientifiques 
^ he Bois-Marie, 35, route de Chartres, 91440 Bures-sur-Yvette, France 

Laboratoire de Physique Theorique de VEcole Normale Superieure, 



24 rue Lhomond, 75231 Paris, France 



4 Centre de recherches mathematiques, Universite de Montreal 
^j- C.P. 6128, Succ. centre-ville Montreal, Que, H3C 3J7, Canada. 

G\ 

(N 

Q Abstract 

In a previous paper, we presented a matrix model reproducing the topological 
string partition function on an arbitrary given toric Calabi-Yau manifold. Here, we 
study the spectral curve of our matrix model and thus derive, upon imposing certain 
minimality assumptions on the spectral curve, the large volume limit of the BKMP 
"remodeling the B-model" conjecture, the claim that Gromov-Witten invariants of 
any toric Calabi-Yau 3-fold coincide with the spectral invariants of its mirror curve. 
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1 Introduction 

In a previous paper pQ, we presented a matrix model that computes the topological string 
partition function at large radius on an arbitrary toric Calabi-Yau manifold X. The goal 
of this paper is to determine the corresponding spectral curve S. 

That the partition function of a matrix model can be recovered to all genus from its 
spectral curve was first demonstrated in [2]. [3] pushed this formalism further, showing 
that symplectic invariants F g (S) can be defined for any analytic affine curve S, with no 
reference to an underlying matrix model. These invariants coincide with the partition 
function of a matrix model when S is chosen as the associated spectral curve. The sym- 
plectic invariants F g satisfy many properties reminiscent of the topological string partition 
function jU [5], [7] , motivating Bouchard, Klemm, Marino, and Pasquetti (BKMP) [8], 
building on work of Marino [9J, to conjecture that F g (S) in fact coincides with the topo- 
logical string partition function on the toric Calabi-Yau manifold with mirror curve S. 
BKMP successfully checked their claim for various examples, at least to low genus. The 
conjecture was subsequently proved in numerous special cases [TO ] fTT | H2 | [13 ] Ibl ] fT5] . 

Bouchard and Marino [16] noticed that an infinite framing limit of the BKMP conjecture 
for the framed vertex, X = C 3 , implies a conjecture for the computation of Hurwitz 
numbers, namely that the Hurwitz numbers of genus g are the symplectic invariants of 
genus g for the Lambert spectral curve e x = y e~ y . This conjecture was proved recently by 
a generalization of [ID] using a matrix model for summing over partitions [17], and also by 
a direct combinatorial method [18]. Matrix models and the BKMP conjecture related to 
toric Calabi-Yau geometries arising from the triangulation of a strip were recently studied 
in [IS]. 

In this paper, we derive the large radius limit of the BKMP conjecture for arbitrary toric 
Calabi-Yau manifolds, but with one caveat: to determine the spectral curve of our matrix 
model, we must make several minimality assumptions along the way. To elevate our 
results to a rigorous proof of the BKMP conjecture, one needs to establish a uniqueness 
result underlying our prescription for finding the spectral curve to justify these minimal 
choices. Such a uniqueness result does not exist to date. 

Recall that in [TJ, we first compute the topological string partition function on a toric 
Calabi-Yau geometry Xo which we refer to as fiducial. We then present a matrix model 
which reproduces this partition function. Flops and limits in the Kahler cone relate Xo 
to an arbitrary toric Calabi-Yau 3-fold. As we can follow the action of these operations 
on the partition function, we thus arrive at a matrix model for the topological string on 
any toric Calabi-Yau 3-folds. Here, we follow the analogous strategy, by first computing 
the spectral curve of the matrix model associated to Xo, and then studying the action of 
flops and limits on this curve. 

The plan of the paper is as follows. In section [2] we introduce the fiducial geometry Xo 
and its mirror. The matrix model reproducing the partition function on Xo, as derived 
in pQ, is a chain of matrices matrix model. It is summarized in section [3] and appendix 
[A] We review general aspects of this class of matrix models and their solutions in section 



|4j In section |5j we determine a spectral curve which satisfies all specifications outlined in 
section [4], and demonstrate that it coincides, up to symplectic transformations, with the 
B-model mirror of the fiducial geometry. While in our experience with simpler models, 
the conditions of section [4] on the spectral curve specify it uniquely, we lack a proof of this 
uniqueness property. We thus provide additional consistency arguments for our proposal 



for the spectral curve in section |5~5j Flops and limits in the Kahler cone relate the fiducial 
to an arbitrary toric Calabi-Yau manifold. Following the action of these operations on 
both sides of the conjecture in section [6] completes the argument yielding the BKMP 
conjecture for arbitrary toric Calabi-Yau manifolds in the large radius limit. We conclude 
by discussing possible avenues along this work can be extended. 



2 The fiducial geometry and its mirror 

2.1 The fiducial geometry 

In PQ, we derived a matrix model reproducing the topological string partition function 
on the toric Calabi-Yau geometry X whose toric fan is depicted in figure [T] We refer to 
Xq as our fiducial geometry; we will obtain the partition function on an arbitrary toric 
Calabi-Yau manifolds by considering flops and limits of Xq. 
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Figure 1: Fiducial geometry Xq with boxes numbered and choice of basis of H2(%o,Z). 

We have indicated a basis of H2(Xq,7j) in figure [l| Applying the labeling scheme in- 
troduced in figure [2j the curve classes of our geometry are expressed in this basis as 
follows, 



r i,j 
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It proves convenient to express these classes as differences of what we will refer to as 
a-parameters [I], defined via 
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Figure 2: Labeling curve classes, and introducing a-parameters. 



2.2 The mirror of the fiducial geometry 



The Hori-Vafa prescription [20J allows us to assign a mirror curve to a toric Calabi-Yau 
manifold. Each torically invariant divisor, corresponding to a 1-cone p G S(l), is mapped 
to a C* variable e~ Yp . These are constrained by the equation 



pes(i) 



0. 



Relations between the 1-cones, as captured by the lattice Ah introduced in section (2.1) 
of [1], map to relations between these variables: for a G S(2), 



J2 a»y p = w a 

pes(i) 



(2.1) 



The W a are complex structure parameters of the mirror geometry, related to the Kahler 
parameters w a = Tij, Sij, . . . introduced in the previous subsection via the mirror map, 
as we will explain in the next subsection. 

The Hori-Vafa prescription gives rise to the following mirror curve Cx of our fiducial 
geometry X , 

n+l m+1 

EE^ =0 ' ( 2 - 2 ) 

i=0 j=0 

We have here labeled the 1-cones by coordinates (i,j), beginning with (0, 0) for the cone 
(0, 0, 1) in the bottom left corner of box (0, 0) as labeled in figure [TJ and introduced the 
notation 



x i,j ~ e 



Eliminating dependent variables by invoking (2.1) yields an equation of the form 

n+l m+1 



2 = j=0 



^i-,2 i-,3 * 



(2.3) 



Here, 



l—i—j i j 
i;3 — •*'l a '2 ) 



where we have defined 



x o — #0,0 > #i — #1,0 j #2 — #0,1 • 



(xq '■ X\ : x-i) define homogeneous coordinates on CP 2 . The form of the equation is 
independent of the choice of triangulation of the toric diagram. What does depend on 
this choice are the coefficients Qj. It is not hard to write these down for the fiducial 
geometry X with the choice of basis for H 2 (X ,Z) indicated in figure [I] Explicitly, the 
relations between the coordinates of the mirror curve (2.2) are 



#i,0 



#i— l,0#i— 1,1 R i _ 



#0, 



#oj-i#i,j-i s 



'3-2 



X 



I -J 



#i-2,l #l,j-2 

Solving in terms of Xo, X\, x 2 yields the coefficients co,o = Cq,i 

"i-l 



#t-l,j#i,j-l Tj-x,,--! 
X i—X,j—1 

= Ci,o = 1, 



(2.4) 



c ifi = exp 



c j = exp 



y^ji - k){R k -i + T fc _i i0 

,fc=i 

'i-i 

^(j-O^i + To,^) 
.1=1 



and for i,j>0 



-»J 



exp 



j-i i-l j-i j-i 

+ j - l)T ,o + $> - fc)(i2 fc _i + T M ) + $>' - 0(^,-1 + T 0>1 ) + J2J2 Tk >< 
fc=i i=i fc=i /=i 



Note that the number of coefficients qj, up to an overall rescaling, is equal to the number 
of independent curve classes r t , Sj, Uj. 

In [21] , the thickening prescription was put forth for determining the genus and number of 
punctures of the mirror curve: one is to thicken the web diagram of the original geometry 
to obtain the Riemann surface of the mirror geometry. The procedure is illustrated in 



figure pi We will now verify this procedure by studying the curve (2.3) explicitly. 



Let's consider the curve (2.3) for a single strip (i.e. n = 0) of length m + 1, 



x, 



m+2 




„m+l 



771+1 



mjl 



„m-XJl, 



_m— 1_3 



+ Xq^Xi + £q £ 2 + Ci,i XqX 1 X 2 + C 2 ,o x q x I + C 2 ,l Xq x\x 2 + C 3 , Xq xl + . 



+ Cn+1,0 x x? +1 + C m+ i,i x? +1 x 2 = . (2.5) 



Note that the equation is of degree m + 2, but the point (0 : : 1) is an m + 1 -tuple 
point. By choosing the coefficients to be generic, we can arrange for this singular point 
to be ordinary. The genus formula then yields 



(d-l)(d-2) m(m + l) 










Figure 3: Example of the thickening prescription: depicted are the fan for 0(— 3) — > P 2 , the corre- 
sponding web diagram, and the mirror curve obtained via the thickening prescription. 



In terms of the physical variables Y$, any point on the curve with a vanishing homogeneous 



coordinate corresponds to a puncture. The punctures on the curve (2.5) thus lie at 

(0: 



(0:1 
:i:x\ 



1) 
0) 
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ra + 1 
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where x\, % 



(1:0: -1) 

,771 + 1, are the solutions of the equation 



1 + X\ 



m 



/ J CtiXi 
0=1 



i+1 







Note that we reproduce the 2m + 4 punctures expected from the thickening prescription 
of the toric diagram. 

For the general case parametrized by (m, n), the degree of the curve is d = m + n + 2, and 
we have an ordinary m + 1-tuple point at (0 : : 1) and an ordinary n + 1-tuple point at 
(0:1:0). The genus formula now yields 

(m + n)(m + n + 1) m(m + 1) n{n + 1) 







The punctures lie at 



mn . 



(0: 
(0:1 

(1 : x\ 



1) 
0) 

0) 



1:0: xi) 



m + 1 

71 + 1 

m + 1 
n + 1, 



with x^ the roots of Yl^Lo c ifi x \ = an d ^2 the roots of Y^j=o c o,j x i = 0- Again, we see 
that we reproduce the thickening prescription. 
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2.3 The mirror map 

Above, we have distinguished between Kahler (A-model) parameters w a and complex 
structure (-B-model) parameters W a . At large radius/complex structure, these are identi- 
fied between mirror pairs, but this identification is corrected by the so-called mirror mapjj 

W a = w a + 0{e- w °). (2.6) 

The exponentials of the parameters W a appear as coefficients in the equation defining the 
mirror curve. They are global coordinates on the complex structure moduli space of the 
mirror curve. To compare expressions obtained in the A-model to those obtained in the 
B-model, all expressions are conventionally expressed in terms of flat coordinates w a . On 
the A-model side, these coordinates enter (in exponentiated form denoted generically as 
Q a ^ below) in the definition of the topological vertex. On the B-model side, they arise 
as the appropriate periods of a meromorphic one-form A, defined in terms of the affine 



variables x = — , y = — in the patch xq ^ of the curve (2.3) as 



i dx 
A = logy—. 

x 
By calculating these periods as a function of the coefficients defining the mirror curve, we 



obtain the mirror map (2.6) 



The coordinates w a are not globally defined functions on the complex structure moduli 
space. In the slightly clearer compact setting, this is due to the fact that the symplectic 
basis {o:a, (3 A } of H 3 (X,Z) in which we expand Q (the compact analogue of the mero- 
morphic 1-form A introduced above) such that the coefficients of a a furnish our (local) 
coordinate system of the complex structure moduli space, undergo monodromy when 
transported around a singularity in moduli spacejj A good choice of coordinates in the 
vicinity of a singular divisor D hence involves a choice of basis forms that are invariant 
under monodromy around that divisor. 



-'^One could take exception to this nomenclature, arguing that the parameters W a are the geometric 
parameters on both sides of the mirror, and refer to the w a as the instanton or quantum corrected 
parameters. In such conventions, the curve classes in the various toric diagrams should be labeled by 
upper case letters. 

2 Note that the symplectic basis makes no reference to complex structure, one might hence be led to 
believe that a global choice (i.e. one valid for any choice of complex structure) should be possible. This 
is not so. We consider the family 7r : X — > S, with S the complex structure moduli space. The fiber 
over each point w £ <S, w w = X w , is the Calabi-Yau manifold with the respective complex structure. 
H n (X w ,C) fit together to form a vector bundle Tq over S, with a canonical flat connection, the Gauss- 
Manin connection. Using this connection, we can parallel transport a symplectic basis of H (X W ,C) 
along a curve in S. As S is not generically simply connected (due to the existence of degeneration points 
of the geometry) , this transport may exhibit monodromy. Note that £1 can be defined as the section of 
a sheaf in the Hodge filtration of H 3 which extends to the singular divisor, hence is single valued. The 
monodromy in our choice of flat coordinates is therefore entirely due to the choice of symplectic basis. 



3 Our matrix model 

We derived a chain of matrices matrix model that reproduces the topological string par- 
tition function on Xq in [1]. For X of size (n + 1) x (m + 1), as depicted in figure [IJ it is 



given by 

m+l „ m+l „ 

ZMu(Q,9s,^m+i,^o) = A(X(a m+1 )) A(X(d )) ][ / dM { TT / dR, 

i=0 JH N (ri) i=1 Jh n (r+) 

m m 

TT e ^ tr [^W-\-i ^ M ')] TT e^ tr [^-iW-J-^W-')] 

8=1 

m+l m 



itr^-Afi.i)^ Y[ e (Si+^)ti bxMi 



i=l t=l 



e trln/o(M ) tr ln/ m+ i(JW m+ i) 



jjetrlnAW). (3.1) 



i=l 



We give the explicit expressions for the various functions entering in this definition in 
appendix |A| Here, we briefly explain some of its general features. 



The matrix model (3.1) is designed to reproduce the topological string partition function 
on the toric Calabi-Yau manifold 3Lq as computed using the topological vertex [22]. Recall 
that in this formalism, the dual web diagram to the toric diagram underlying the geometry 
is decomposed into trivalent vertices. Each such vertex contributes a factor C(cti, otj, «&) 
[22], where aij denote Young tableaux (partitions) of arbitrary size, one associated to each 
leg of the vertex. Legs of different vertices are glued by matching these Young tableaux 
and summing over them with appropriate weight. 

Aside from the coupling constant q s and Kahler parameters of the geometry, denoted 
collectively as Q, the matrix model (J3JJ) depends on partitions c?o, a m +i associated to 
the outer legs of the web diagram, which we choose to be trivial in this paper. The two 
classes of integrals dRi and dMi correspond to the two steps in which the topological string 
partition function on the fiducial geometry £ can be evaluated: First, the geometry can 
be decomposed into m + l horizontal strips, with partitions Oj,j+i and <x, : i associated 
to the upper and lower outer legs of the associated strip web diagram, j = 0, ...,n 
counts the boxes in figure [T] in the horizontal direction, i = 0, . . . , m + 1 is essentially 
the strip index. Each such strip has a dRi integration associated to it. The partition 
function on such strips was calculated in [23]. Following |12j . we introduce two matrices 
Mi, Mj+i per strip. Their eigenvalues encode the partitions ct^ and <x,-,j+i for all j. To 
work with finite size matrices, we introduce a cut-off d on the number of rows of the 



Young tableaux we sum over. As we argue in section 5.2.1 , our matrix model depends on 
d only non-perturbatively. The strip partition function is essentially given by the Cauchy 
determinant of the two matrices Mi, M i+1 pQ, and the dRi integrals are the associated 
Laplace transforms. Gluing the strips together involves summing over the partitions o^. 
This step is implemented by the dMi integrations. To obtain a discrete sum over partitions 



from integration, we introduce functions fi(Mj) with integrally spaced poles. Integrating 
Mi along appropriate contours then yields the sum over partitions as a sum over residues, 
the potentials Vs t chosen to provide the proper weight per partition. 

4 Generalities on solving matrix models 

4.1 Introduction to the topological expansion of chain of matri- 
ces 

Chain of matrices matrix models have been extensively studied (see Mehta's book [21] and 
the review article [22]), and the computation of their topological expansion was performed 
recently in [261 1211 • 

The solution provided in [27] is based on the computation of the spectral curve 5mm of 
the matrix model. In J27J [26], only the case of potentials whose derivatives are rational 
functions is considered, and similarly to the one matrix model, the planai|j expectation 
value of the resolvent of the first matrix of the chain is shown to satisfy an algebraic 
equation. The spectral curve is defined to be the solution locus of this equation. A general 
recipe is provided in [261 EZ] to obtain the spectral curve from algebraic equations and 
analyticity properties related to rational potentials and integration contours. Here, our 
potentials contain logs of ^-functions. As they are not rational, we will have to present 



a slight extension of the recipe of [27] in section 4.2.2 This extension from rational 
potentials to analytical potentials, although not published, is straightforward, and the 
derivation of these results will appear in [28]. In some sense, the derivative of \ag(x) 
can be viewed as a rational function with an infinite number of simple poles, i.e. as a 
limit of a rational function. More precisely, as an expansion in powers of q, to each order, 
it is a rational function. Since the spectral curve can be described by local properties, 
independent of the number of poles, one can take the limit of the recipe of [261 127]. This 
is what we shall do in section 14.2.21 below. 

Having found the spectral curve 5mm of the matrix model, we will compute its symplectic 
invariants 

F g (SuM) , 9 — 0, 1,2,3, .. . 

Symplectic invariants F g (S) can be computed for any analytical plane curve S, and thus 
in particular for S = Smm- For a general S they were first introduced in |3], as a 
generalization of the solution of matrix models loop equations of [2]. Their definition 



3 For matrix models with iV-indcpendent polynomial potentials whose g s dependence is given by an 
overall prefactor, the planar limit coincides with the large N limit, but this correspondence can fail if the 
potential or the integration contours have a non-trivial N or g s dependence. The planar limit is defined 
by keeping only planar graphs in the Feynman graph perturbative expansion around an extremum of the 
potential. However, it is helpful to have in mind the intuitive picture that the planar limit is similar to 
a large N limit. 



is algebraic and involves computation of residues at branch points of S. We recall the 
definition below in section I4.3L 



4.2 Definition of the general chain of matrices 

We consider chain of matrices matrix models of the form 



Z= fdM 1 ...dM L e-^ ^ «<*» e^ E - C > M < M >^ 



(4.1) 



Note that aside from the potentials Vi(Mi), the only interactions are between nearest 
neighbors, whence the name "chain of matrices." Chain of matrices matrix models can be 
solved when the interaction terms between different matrices are of the form Tr MiM i+1 , 
as is the case here. 

£ can be any ensemble of L normal matrices of size N x N, i.e. a submanifold of C LN 
of real dimension LN 2 , such that the integral is convergent. £ can be many things; for a 
chain of matrices model, it is characterized by the contours on which eigenvalues of the 



various normal matrices are integrated (see [29] for the 2- matrix model case). For (4.1 ) to 
have a topological expansion, £ must be a so-called steepest descent ensemble (see [SO] . 
section 5.5). For a generic ensemble £ which would not be steepest descent, \nZ would 
be an oscillating function of l/g s , and no small g s expansion would exist, see [3l] . 

The matrix model introduced in pQ and reproduced in section [3] was defined to reproduce 
the topological string partition function, which is defined as a formal series in g s , and 
therefore has a topological expansion by construction. 

An ensemble £ is characterized by filling fractions rij^, 

L 

£ = \{£, , Si = Hsfrtf x 7 £' x • • • x 7 2i'«) , (4.2) 

4 = 1 

where H^iji 1 x ■ ■ ■ x 7 £*) is the set of normal matrices with rt\ eigenvalues on path 7 i, 
rii eigenvalues on path 72 , . . . , nk eigenvalues on path 7 £. 

As the filling fractions rij^ must satisfy the relation 



E 



n j,i = N 



for all i, only Y2i(^i — 1) of them are independent. 

We also allow some paths 7^ to have endpoints where e^Sid^W)" 1 '^ 1 ) ^ 
indeed, in our matrix model, the matrices Ri are integrated on H(M.+ ). 
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4.2.1 The resolvent 

The spectral curve encodes all Wi(x), the planar limits (see footnote |3| of the resolvents 
of the matrices Mi, 

Wi{x) = g s ( tr 



X ^i I planar 



see equation (4.5) below. The respective Wi can be expressed as the Stieljes trans- 
form 

= f PiWW 

y ' J x-x' 
of the planar expectation value of the eigenvalue density Pi(x) of the matrix Mi, 

Pi(x) =g s (tr5(x-Mi)) planar . 

By general properties of Stieljes transforms, singularities of W{(x) coincide with the sup- 
port of the distribution pi(x)dx: 

• Simple poles of Wi(x) correspond to delta distributions i.e. isolated eigenvalues. 

• Multiple poles correspond to higher derivatives of delta distributions. 

• Cuts correspond to finite densities, the density being the discontinuity of Wi(x) 
along the cut, 

Pi{x) = ^- {Wi{x - zO) - Wi{x + zO)) . (4.3) 

In particular, cuts emerging from algebraic singularities (generically square root singular- 
ities) correspond to densities vanishing algebraically (generically as square roots) at the 
endpoints of the cut. Cuts emerging from logarithmic singularities correspond to constant 
densities. 

4.2.2 The spectral curve of the general chain of matrices 

When all V( are rational, the spectral curve was found in [261 [27], an d it is algebraic. We 
present here a generalization of this result to more general potentials. The derivations of 
these results will appear in [28]. 

The spectral curve can be obtained by the following procedure: 
1. Consider a compact Riemann surface C of genus 

L 
= 5>,-l), 
t=l 



where ki denotes the number of cuts of the z-th matrix, as implicitly defined in (4.2 ) 
2. Look for L + 2 functions on C, 

Xq(z), xi(z), x 2 (z), ..., x L (z), x L+1 (z) : C - 

11 



The Xi are to be holomorphic away from points z G C at which V^ixi-i) or 
Vl +1 (x.i + i) become singular, and satisfy the functional relations 

c i _ 1 x i _i(z) + CiX i+1 (z) = V-(xi(z)) . (4.4) 



Recall that the q are the coefficients of the interaction potentials in (4.1 ). We have 

set Co = cl = 1. 

For each i = 1, . . . , L, the Riemann surface C can be realized as a branched covering 

of CP 1 by the projection xi : C — > CP 1 . A choice of branched covering is not unique: 

the choice consists in the set of cuts connecting branch points (recall that these are 

points at which dxi(z) = 0). We will determine an appropriate covering below in 

step[3J 

3. If some path 7^ has an endpoint a (called "hard edge" in the matrix model litera- 
ture, see [32J), then choose a pre-image Oj G x~ 1 {a) and require 

dxi{a,i) = and Xi-i(z) has a simple pole at z = ai. 

The topological recursion is proved in [27] without hard edges, but it is not difficult 
to see, by mixing the results of [32], [33] and [27], that the topological recursion con- 
tinues to hold in the presence of hard edges. The proof will appear in a forthcoming 
publication [2H]- Here, we shall assume that it holds. 

4. Choose some contours Ajj, j = 1, . . . , fc« in CP 1 , such that each Aj t i surrounds all 



points of the contour 7^ (related to the matrix ensemble £{ defined in (4.2)) in the 
clockwise direction and no other contour 7^. For x G CP 1 not enclosed in the 
contours Aj t i, j = 1, . . . , k iy and given a connected component Aj,% of the pre-image 



of the contour Aj,i under Xi 



ii* i \"^j,i 



define the function 



lK ! 2i7i 4-f J a.. x-xAz) y ! 



i=i 



A 



Generalizing [26] to non-polynomial potentials, we claim that a choice of Aj,i exists 
such that Wi(x) is the planar limit of the resolvent of the matrix Mj. In the following, 
it is this choice that will be referred to as A^%. 

Notice that not all Ajj will be homologically independent on C. We require that 
we have $j = J2i=i(^i ~ 1) homologically independent Aj/s, which coincides with 
the genus of C. As a condition on the choice of branched covering, we impose that 
Aj } i and Oj lie on the same sheet of Xj. This condition, in our experience, uniquely 
fixes this choice. We will assume that this is the case. We refer to the sheet of Xi 
containing Aj,i and Oj as the physical sheet for x%. 
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5. In accord with (4.3), we consider the discontinuity of Wi(x) along the j-th cut. It 
is given by 



Disc Wi(x) 

j 



as we explain in figure [5j 



— (Wi{x + ) - Wi{xJ)) 

-— Ci-i Disc Xi-i , 
2m j 



(4.6) 




Figure 4: The integration contour Aji on the x-plane, and its image Aji on C. 









Figure 5: The preimage of the points x+ and X- of (4.6| are depicted as dots in the above diagram, Aj t 



■*3>« 



is given by the blue contour, and the preimage of the cut is drawn in red. To take the limit x+ — > X—, 
one must first shift the contours. The second and fourth term on the RHS of the above diagrammatic 



equation then cancel, yielding the RHS of (4.6 1 



The definition (4.2) of the matrix ensemble E{ is the condition that there are n 



j-i 



eigenvalues of M, on the contour 7^, hence corresponds to imposing the filling 
fraction conditions 



for i — 1, 



L, j = 1, 



2ixi 



, fcj. 



A, 



In our experience, the conditions enumerated above have a unique solution and define a 
unique spectral curve. As emphasized in the introduction, a formal uniqueness proof is 
however still lacking. 

The spectral curve is defined as the data of the Riemann surface C, and the two functions 
Xx(z) and £2(2), 



Smm = (C,x 1 ,x 2 ). 



(4.7) 
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4.3 Symplectic invariants of a spectral curve 

Once we have found the spectral curve Smm of our matrix model, we can compute the 
coefficients F g in the topological expansion of its partition function, 

oo 

in z = Y,g 2 r 2 F 9 , 

9=0 

by computing the symplectic invariants of this curve, 

Fg — F 9 (Smm) , 
following [27] . 
Let us recall the definition of these invariants for an arbitrary spectral curve S. 

Let S = (C, x, y) be a spectral curve, comprised of the data of a Riemann surface C and 
two functions x(z), y(z) : C — > C, meromorphic on C away from a finite set of points 
(we wish to allow logarithms) n We will assume that dx is a meromorphic form on all of 
C. 

4.3.1 Branchpoints 

Let cij be the branch points of the function x, 

dx(di) = 0. 

We assume that all branch points are simple, i.e. that dx has a simple zero at aj. This 
implies that in the vicinity of a,, the map x is 2 : 1. We introduce the notation z ^ z 
such that 

x(z) = x(z). 

z is called the conjugate point to z, and it is defined only in the vicinity of branch points, 
as depicted in figure [6j 

We also require that the branch points of x and y do not coincide, such that dy(a,i) ^ 
and y{z) therefore has a square-root branchcut as a function of x at x(aj). If y is finite 
at <2j, its local behavior is hence given by 

y(z) ~y(di) + Ci^Jx(z) -x{di) . 

If di corresponds to a hard edge, we require y to have a pole here. Its local behavior is 
hence given by 

y(z) ~ -. . 

^x(z) - X{di) 



4 In fact, the most general setting in which this formalism is valid has not yet been established. We 
will state it within the generality we need here i.e. we assume that dx is meromorphic forms on C (this 
allows x and y to have logarithms). 
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Figure 6: At a regular branch point a G C of x, y as a function of x has a branchcut y ~ y(a) + 
C\Jx — x(a). If z is a point on one branch near a, we call z the conjugate point on the other branch; it 
has the same x projection, x(z) = x(z). Notice that z is defined only locally near branch points. If we 
follow z from ai to a 2 , z may have to jump from one branch to another. 

4.3.2 Bergman kernel 

On a curve C, there exists a unique symmetric 2-form B(zi, z 2 ) with a double pole on the 
diagonal z\ = z 2 and no other poles, with the following normalization on ^4-cycles, 

B( Zl ,z 2 ) = 0. 
In any local coordinate near Z\ — Z2, one has 



B(z 1 ,z 2 ) ~ 



dz\ dz2 



+ regular . 



y Zi - z 2 ) 
B is called the Bergman kernel of C, or the fundamental 2-form of the second kind 



4.3.3 Recursion kernel 



We now define the recursion kernel K as 



K(z ,z) 



// B(zo, z>) 



2(y(z)-y(z))dx(z)- 



This kernel is a globally defined 1-form in the variable zq e C. In the variable z, it is the 
inverse of a 1-form (that means we have to multiply it with a quadratic differential before 
computing any integral with it); it is defined only locally near branch points of x, such 
that K(zq,z) = K(zq,z). At the branch points, it has simple poles, 



K(z ,z) 



B(zq,z) 



2dx(z) dy(z 



+ regular. 
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4.3.4 Topological recursion 

Correlation forms W n (z\, . . . ,z n ) (not to be confused with the resolvents Wi(z) intro- 
duced above) are symmetric n-forms defined by 

Wi°\z) = -y(z)dx(z), 

W^°\z 1 ,z 2 )=B(z 1 ,z 2 ), 
and then by recursion (we write collectively J = {z\, . . . , z n }), 

i 

+ EE<|(^)<- ) |/|(^\/) 
h=0 ICJ 

where Yli * s the sum over all subsets of J, restricted to (h, I) ^ (0, 0) and (h, I) ^ 
(9, J)- 

Although it is not obvious from the definition, the forms W n g are symmetric. For 2 — 
2g — n < 0, they are meromorphic n-forms with poles only at branch points. These poles 
are of degree at most 6g — 4 + 2n, and have vanishing residues. 

For the one matrix model, the Wn coincide with the n-point function of the trace of the 
resolvent at order g in the topological expansion. 

4.3.5 Symplectic invariants 

Finally, for g > 2, we define the symplectic invariants F g (also denoted Wo(g) in [3]) 
by 

F 9 (5) = -^-^Res$(.)< ) (^) ) 

3 i 

where $ is any function defined locally near branch points of x such that <i$ = ydx. 

The definitions of F\ and F are more involved and we refer the reader to [3] . F Q is called 
the prepotential, and i*\ is closely related to the determinant of the Laplacian on C with 
metrics \ydx\ 2 , see [351 l36j . 

The F g (S)'s depend only on the orbit of S under the group of transformations generated 
by 
£H : S i-)- S = (C, x, y + R(x)) where R(x) is any rational function of x, 

5 : S i— > S = (C, f(x),y/f'(x)) where f(x) is an analytical function of x, with /' 

rational, such that df = f'dx has the same number of zeroes as dx, 

6 : S h-» S = {C,y, -x). 
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These transformations are symplectic, i.e. they leave dx A dy invariant. 
The symplectic invariants are homogeneous of degree 2 — 2g, 

F g (C,x,Xy) = X 2 - 29 F g (C,x,y). (4.8) 

In particular, they are invariant under the parity transformation F g (C, x, —y) = F g (C, x, y). 

5 The spectral curve for the topological string's ma- 
trix model 



Applying the procedure outlined in section 4.2.2 to our matrix model, we will determine 



a spectral curve Smm(%o) in this section. [27] demonstrated that for a chain of matrices, 
we have _ 

\nZ = J29 2 r 2 F g (S MM ), 

9 

with F g the symplectic invariants of [3j. In our case, since we have engineered our matrix 
model to yielcjj GW g (Xo) as its partition function, re-computing the partition function 
via the methods of [27] will yield 

GW g (X Q ) = F g (S MM ) ■ 

This relation is already quite interesting, as it allows for explicit computation of the 
Gromov-Witten invariants. Our goal however will be to go further. We will argue that 
iS>mm is symplectically equivalent to the mirror spectral curve S% o of section 



2.2 



Since the F 5 's are symplectic invariants, this will imply the BKMP conjecture for jE > 
i.e. 

GW g (X ) = F g (S Xo ). 

5.1 Applying the chain of matrices rules 



We now apply the rules of section |4.2.2| to the chain of matrices model introduced in 
section HI 

• Recall that the integration ensembles for the matrices Mq and M m+ \ are such that 



for each matrix, all eigenvalues are integrated on the same contour (A.7). Hence, 
ko — k m+ i = 1, and the corresponding filling fractions are equal to N. For i = 
1, ... ,771, the matrix Mi is integrated on Hipf^ x 7^ x • • • X 7^), where 7^ is a 



5 As we have here reserved the notation F g for the symplectic invariants of our matrix model, we refer 



to the topological string free energies as GW g . 
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contour which surrounds all points of the form g a J. i+N . There are thus ki = n + 1 
filling fractions, each equal to d. The matrices Ri are integrated on H(R~[). We 
denote the number of their cuts by ki. Hence, ki = 1, with the respective filling 
fraction equal to N. 

According to condition [I] of section 4.2.2 , the genus of the spectral curve C is thus 
given by 



m+l 



m+l 



Q = j2(k l -i) + J2^- 1 ) 



nm . 



i=0 



t=l 



Following conditionp^of section 4.2.2, we introduce functions Xi(z), i = 0, . . . , m+l, 
associated to the matrices Mi, and functions yi(z), i — 1, . . . , m + 1, associated to 
the matrices Ri, as well as two additional functions yo{z) and y m+ 2(z) at the ends 
of the chain. 

They must satisfy the following requirements: 



Since there is no potential for the matrices Ri, equation (4.4) implies that we 
have, for i — 1, . . . , m + 1, 

Xi(z) - Xi_i(z) = 0. 

We can hence suppress the index i on these functions, x(z) = Xi(z). 
For i — 1, 



,m, equation (4.4) gives 



yi(z)-y i+ i(z) = 2Vi(x(z))-Vi i+1 (x(z))-Vi i _ 1 (x(z))-g 
and 



fii x ( z )) g s Si + iii 



y (z) - yi(z) = Va o (x(z)) - V^(x(z)) - g t 



y m +i(z) - y m+ 2{z) = Vi m+l {x{z)) - Vi m {x(z)) - g. 

More explicitly, in terms of the function 

i/j q (x) = xg'(x)/g(x) , 
whose small g s expansion 



fi(x(z)) x(z) 

f'Mz)) 
fo(x(z)) ' 

fm+l(%( Z )) 



f m+ i{x{z)) ' 



i) q {x) 



1 oo 

--Y 

lng z — ' 



-l) n B„ 



n=0 



IV. 



(\nq) n Lh.^l/x) 



no L 



, >(1-- 
mg L x' 



lng 



2(x - 1) ^^ (2n)\ 

v ' n=l v ' 



E 



B 



2n 



(InqyoLi^ix) 



we worked out in appendix A of [T], we obtain 

x(z)(y i+1 (z) -Vi{z)) 
= i-K + g s Si - g s YtWiWM*)) ~ ^ q (q aj ' i+1 /x(z)) - i> q {q a ^/x{z))) 



■j 



+9s fM*)) ' ( J 

as well as 

x{z)( Vl (z)-y (z)) = -g. y %2if> q {<f'' /x{z))+g s y %2il> q {q a * 1 /x{z)) 



x(z) 



3 

d-1 

^—^ ^-— ' x(z) — q a i,o+k ' 

x(z)(y m+2 (z) - y m+1 {z)) = -g s J2Ml aj ' m+1 / x ( z )) + 9s ^2^ q {f s ' m /x{z)) 

j j 

d-i , , 

-q W X ^' 

j k=o v ' H 

Note that we have explicitly used the fact that the partitions <x, jm+ i and o^o 
are chosen to be trivial. 

Since the integral over Ri is over H N {R + ), i.e. its eigenvalues are integrated on 
R+, the integration contour has an endpoint (hard edge) at yt = 0. Condition 
Whence requires that at a pre-image y i ~ 1 (0), which we will refer to as ooj, the 
following holds 

Vi(oOi) — , dyi(ooi) = , x(z) has a simple pole at z = ooj . 

Furthermore, introducing a local parameter 2 in the neighborhood of oo», the 
above translates into 

yi(z) ~ z 2 , x(z)~l/z. 

Hence, Vi = 1, . . . , m + 1, 

y t ~ O(l/o; 2 ) . 



The relations (5.1) imply that near ooj, we have 



x(yj+i-Vj) ~ i-K + g s Sj + g s 5^(2a u - ojj+i - ojj_i) + 0(l/x). 



«=o 



In particular, it follows that ooj 7^ ooj. Thus, all points {001, . . . , oo m+ i} C x 1 (oo) 
are distinct , i.e. condition M requires that 2 _1 (oo) have at least m + 1 points. 
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m+1 






^ ctij) 4xb 



<5^> <2^ 






Figure 7: The spectral curve of our matrix model can be represented as follows. The cover of CP 1 
provided by x has to + 1 sheets. Instead of the projective plane of x, we represent the sheets of In a;, 
which are cylinders. Cycles A.^% appear in sheets i — 1 and i. They enclose singularities of the resolvent 
Wi. Algebraic cuts are represented as vertical cylinders, and poles and log singularities are represented 
as grey strips. There is only one cycle Aj t o (which is in sheet 0) and one Aj >m +i (in sheet to), and they 
enclose only poles or log singularities of j/o resp. y m+ i. 



We will make the minimal assumption that x _1 (oo) has exactly m + 1 elements 
that are simple poles of x, and that x has no further singularities, i.e. that x is a 
meromorphic function of degree m + 1 on C. 

By condition [5j since for i — 0, . . . ,m + 2 there are d eigenvalues of Mi of the form 
g a j,i+ N surrounded by the path Ajj , we have the (m + 2) x (n + 1) filling fraction 
conditions 



1 

2iTT 



Hi dx = d g s for i — 0, . . . , m + 1, j — 0, 



n. 



A, 



x hence defines anm + 1 sheeted cover of CP 1 . Considering the function In re instead, 
with singularities at x = and x = oo, each sheet of this cover is mapped to a cylinder. 
We have depicted this covering in figure [7j and indicated the singularities of yi on each 
sheet: algebraic cuts are represented by vertical cylinders, and poles and logarithmic cuts 
by grey strips. 

In sheet i we have represented some contours Aj,i whose image under the projection 
x : C — > CP 1 surrounds all points of type g a -? i+N . 

For i = 1, . . . ,m, the resolvent Wi(x) of the i th matrix Mi is computed as a contour 
integral around the sum over j of cycles A.j,% on sheet i, 



Wi(x) 



n 

y — 

^ 2in 



j=0 



yi(z')dx(z') 



A, 



X 



X(Z' 
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Also, as argued in [T] , the potentials of M and M m+1 are such that in fact the matrices 
M and M m+1 are frozen, and thus their resolvents contain only poles. In terms of the 
functions yo and y m +i, we conclude that the singularities of yo in Aj$ in sheet 1 and the 
singularities of y m +i in Aj, m +i in sheet m + 1 can be only poles, not cuts. 
Since condition [T] requires that the genus be g = nm, we see that there can be no other 
cuts than the ones already discussed - the genus would be higher, otherwise. 

5.2 Symplectic change of functions 

The spectral curve of the matrix model is Smm = (C, x, yo), and our goal is to relate it to 



the mirror curve described in section 2.2 The mirror curve is described via the algebraic 
equation (2.3) in the two functions X\, x 2 '■ C mirror — > CP 1 (in the patch x = 1). We wish 
to obtain a similar algebraic description of C. Due to log singularities in y , to be traced 
to the small g s behavior of ip q {x), an algebraic equation in the variables {x,yo) cannot 
exist (recall that x is meromorphic). In this section, we shall, via a series of symplectic 



transformations on the yi of the type enumerated in section 4.3.5 arrive at functions Y^ 
that are meromorphic on C, and hence each present a viable candidate to pair with x to 
yield an algebraic equation for C. 

Essentially, we wish to introduce the exponentials of yi. While this will eliminate the log 
singularities, poles in yi would be elevated to essential singularities. We hence first turn 
to the question of eliminating these poles. 

5.2.1 The arctic circle property 

On the physical sheet, the interpretation of a pole of yi is as an eigenvalue of the matrix Mj 
with delta function support. Such a so-called frozen eigenvalue can arise in the following 

way: 

The sum over all partitions is dominated by partitions close to a typical equilibrium 
partition, i.e. a saddle point. The typical partition has a certain typical length referred to 
as its equilibrium length n. All partitions with a length very different from the equilibrium 
length contribute only in an exponentially small way (and thus non-perturbatively) to the 
full partition function. Introducing a cutoff on the length of partitions which is larger 
than the equilibrium length hence does not change the perturbative part of the partition 
function. Now recall that when we defined the h^) of a representation 7 in appendix |Al 
we introduced an arbitrary maximal length d such that ^(7) < d and set 

^(7) = a 7 + d-i + -fi. 

Setting 7j = for d > i > n yields hi that do not depend on the integration variables, 
hence are frozen at fixed values. This behavior is referred to as the arctic circle prop- 
erty [37J, as all eigenvalues beyond the arctic circle situated at equilibrium length n are 
frozen. 
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Returning to our matrix model, the eigenvalues of the matrices M t are given by q 



( h 3,i) 



j,i)l 



For d > I > rijj, they are frozen, and thus contribute poles to y t by (4.5) (recall that poles 
of the resolvent correspond to eigenvalues with delta function support) in the physical 
sheet. We will assume that these are the only poles in the physical sheet and we subtract 
them to obtain new functions ft, 



y (z) =^)!/o(2)-VV 



g s x{z) 



^■— ' *—* x(z) — q a 3' Q+k 



d-l , V 

g s x{z) 



y m+2 (z) = x(z)y m+2 (z) + E E x{z) _ qaj , m+1+ k 



and for i — 1, . . . , m + 1 



, • • • , "« i j-, 



~ ( \ _ t \ ( \ _ V^ V^ 9sx{z) ^ ^ g s x{z) 

yi{z) - x{z)yi{z) ^ Z^ x ( z \ _ q a j:i +k + Z^ Z^ x f z \ _ q aj,i-i+k ■ 

We have set 

Tljfl = , % m+ i = . 

Notice that at large x(z) in sheet i we have 

y ~ 0(\/x(z)) , ftn+2 ~ 0(l/x(z)) 
and for z = 1, . . . , m + 1 

ft ~ 9s ^2(nj,i - n jtt -i) + 0(l/x(z)). 



As a general property of ip q , we have for any integer rij^ < d 

d—n,j t i—l 

ipqiq^/x) = ip q (q aj ' i+d - n ^/x) + J2 



x 



k=0 * 

Hence, the loop equations for the new functions ft read 

ft+i(z) -Vi(z) = in+g s Si 

+g s J2(2ip q (q a ^ i+d - n ^/x(z)) - ^ q (q a ^+ 1+d - n ^ l+1 /x(z)) 
j 

-V^- 1+d -^-7^))) + 9s X{ fl^ ] > 

!/!(*)-#>(*) = ^X)V'g(9 aA0+<l AW)-^Z)^(9 ai,1+<l " Bi,1 AW)» 
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a+d-n 




Jti , with rij.i chosen such that 



Figure 8: We shift the cut-off d on the representation lengths, d —>■ n 

frozen eigenvalues in the expected distribution of the hi are suppressed. In the limit of vanishing spacing 

(g a — > 0), the equidistant frozen eigenvalues give rise to a constant eigenvalue density region. 



y m +2{z) - y m+ i(z) = g s J2^i a] ' m+1+d / x ^ ~ 9sJ2^ aj 



l +d—nj 



>/x(z)). 



The nj t i in the above definitions are defined as the equilibrium lengths, i.e. by the property 
that the functions y, have no poles on their physical sheet. That such a choice of n^i exists 
is suggested by the arctic circle property. 

Note that the rij ti can also be specified by the fact that q a ^ i+d ~ n J' i be the beginning of the 
cut encircled by 7^. As we have identified the discontinuities of y^ to lie across branchcuts 
of x, this implies that x has ramification points at the element of x -1 (g aj '' i+ d ~ n i> i ) lying on 
the physical sheets of y{. 

Note that the arctic circle property also implies the perturbative independence of our 
expressions from the arbitrary cut-off d. Changing d to d + d' merely introduces d! new 
frozen eigenvalues h{. This independence from d is important in establishing the equality 



between the topological string partition function and our matrix integral (3.1), as the 
topological vertex formulae in fact are formulated in the limit d — ¥ 00. 



5.2.2 Obtaining globally meromorphic functions 

We have arrived at functions y~i that have no poles on their physical sheet, and are thus 
safely exponentiated there. We wish now to use the loop equations to obtain functions 
which are globally well-behaved. 

To this end, we note that since the Gromov-Witten invariants are defined as a formal 
power series in g s , we can compute the spectral curve order by order in g s , invoking the 
following small In q expansion pQ : 



^ q (q a ^ +d - n ^/x) 



1 



ln(l 



x 



X 



B2n9. 



_. 00 

^5 <2 " )! 



fy^Jji' 3>% 



2(x 



a a 3,i 



Lii_ 2n (q' 



d - n i,i+ a i,i 



/x). 
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The functions f-ffi are completely non-perturbative; one can easily check with the above 
expansion that they can be replaced by to every order in g s . 

Introducing new functions X(z) and Yi by the formulae 

x(z) = q d X(z), 

y (z) = \nY (z), 
y m +2{z) =\nY m+2 (z), 
and for % — 1, . . . , m + 1 

ft« = -* W + E ^F £■■-»-) + ^ S W Li -fa--°"/ x W) 

yields loop equations that are algebraic on their right hand side, 

Y i _ TT (X — q a ^ i+1 ~ n ^ i+1 )(X — q a ^ i - 1 ~ n ^ i - 1 ) 



n 



Y i+1 J-. 1 (X - q a J,i- n J,i) 2 

q 



TT rr 2 ( a j,i-«j,i)-( a i,i + l~™3,i + l)-( a J,i-l~ n j,i-l) 



J 



Y i-r (X - <f ■ 



i,i- n j,i 



(Y-nai,o\ 11 



-.aj t o— (aj t i— n,j : i) 



Fl A 4 (X - g°i.o) " ' ' 



j 



i.e. 



r m+2 - 11 (x-^.^0 I/ 9 ' (5 ' 2) 



V. _ _ Y" /r a J,i-l~ n i,i-l 

11 __ P gs(s 1 +-+Si- 1 ) n„(oj,i-' l i,i)-( |I i,i-i-' l i,i-i) IT - 

Yo lig 1.1 X-g' 



^ 3 ! * ^ J ! * 



J 



Since we have argued that the Yi are holomorphic on their physical sheet, and the ratio 
Yi/Y i+ i is purely algebraic, we conclude that the Yi are meromorphic functions on all of 
C. This was the goal we had set out to achieve. 

Note that the above changes of variables have modified the asymptotics at infinity and 
the integrals over the *4-cycles. More precisely, we have 

Vie[l,m + 1]: \nYi~yi~ g s y^(n j:i -n j: i-i) + 0(-), (5.3) 

ODi ODi *■ ' X 

j 

lnY = y ~O(-), 

oo X 
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kiY m+2 = y m+2 ~0{-\ 

oo X 



The filling fraction equation reads 

1 f Vi(z 



2m J Aj . x{z) 



dx(z) = g s (d - (d - n 3)i )) = g s n j}i . (5.4) 



In terms of Y , these conditions can be rewritten as 

n 1 

lnF ~ -g s {Si H h Si-i) + g s )(aj,i - Oj>-i) + 0{- 

ooi * — ' ' X 

and 
1 f , „ dx l f , „ dx l f , ,r l, ( Y i 



, „dX 1 

lnr, — + — 



1m a X 2m 



2m J a X 2m J A . . \ ^ X - cpu-**,* X - g°i,i-i-»w-i 



i r . „dx 

hiY* — 7- + g s [aj,i -n 






1 /" , „dX 1 /■ , „ /^ dX dX 



*j,» "-j,i/ 



2m J A .. X 

= g s n Jti + g s (a jti - n jyi ) 

5.3 Recovering the mirror curve 

We have argued above that X and Y{, and hence in particular Yq, are meromorphic 
functions on C. There must hence exist a polynomial H(X, Y) such that (see e.g. Theorem 
5.8.1 in |38]) 

H(X,Y o ) = 0. 

The facts that X provides anm + 1 sheeted cover of CP 1 and that Yq may have n + 1 poles 
in its physical sheet imply that the polynomial H has degrees at least (n + 1, m + 1). As 
above, we shall choose the minimal hypothesis that it has exactly these degrees. Thus, 

m+1 n+1 

H(X,Y) = J2Y, H ^ XJYi - ( 5 ' 5 ) 



As we saw in section 2.2 projectivizing a generic polynomial of these degrees (yielding 
a homogeneous polynomial of degree m + n + 2) indeed gives rise to a curve of genus 
q = nm. 

We now need to determine the (n + 2) (m + 2) — 1 unknown coefficients of H (H is defined 
up to a global multiplicative constant). 
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The cycle integrals 

dX 



lnYo -^r = 2i7r^ s o^ 

provide (n + l)m constraints on the coefficients of H. We also have m + 1 constraints for 
the behavior at ooj, i — 1, . . . , m + 1, 

Res In r "v = S'sC'S'i H ^ «Si_i) - # s /",(%,« _ a i,i-0 • 

i=o 

Finally, requiring that Y has poles at q aj -° and 1^+2 has zeroes at q a i> m + 1 gives another 
2(n + 1) constraints, which we may write as 

Res InA — — = o s o,0) 

a J -° ' o 

Res In A — = g s a j)m+1 . 

q a J,m+l Y m+2 

This gives enough equations to completely determine H. Knowing H, we know the 
location of branch points as functions of cij/s and Si's, and can hence determine the rijj 
by requiring that q a i> i ~ n i> i be a branch point. 

Notice that we can choose to express the period integrals in any linear combination of 
^4-cycles. In particular, 

r jv 

(b \ n Y — = 2mg s (a j>i+1 - a jA ) = 1m g s t jti , 

JAj,i+i-Aj,i A 

/ dX 

<b \nY — = 2m g s {a j)i+1 - a j+M ) = 2m g s r j}i . 

JAj^+i-Aj+i,, A 

Similarly, we may also take linear combinations of ^.-cycles together with circles sur- 
rounding the poles or zeroes of x in order to get the Sj t i classes. We hence conclude that 
the periods of the curve H(X, Y ) = yield the quantum corrected Kahler parameters of 
the fiducial toric geometry X Q , allowing us to identify it with the corresponding mirror 
curve. 



5.4 Topological expansion and symplectic invariants 

Following [27], we obtained 

Sum = (C, x, y ) 

as the spectral curve of our matrix model at the end of section 5.1| 



As reviewed in section |4.3[ we can compute the corresponding symplectic invariants 
F 9 (Smm), which assemble to yield the matrix model partition function 



\nZ = Y,9T 2 F g {S MM ) 
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Figure 9: The spectral curve (X, -4= In (lo)) has the following structure: X{z) is a meromorphic function 
of degree m + 1 on a curve of genus q = nm. Therefore it has m + 1 poles and m+1 zeroes. It provides 
a branched covering of CP 1 . We prefer to represent lnX instead of X, and thus we have m + 1 copies 
of the In X-cylinder. In each sheet there is one zero and one pole of X. Yq is a meromorphic function of 
degree n + 1, so that it has n + 1 zeroes in sheet 0, and n + 1 poles in sheet m + 1. We recognize the 
mirror curve 5^ , which is a thickening of the toric web diagram. 



The symplectic transformation ^K of section 4.3.5 maps (C,x,yo) to (C,x, -\nY ) order 
by order in g s . $ maps this to {C,X, j^ lnY ), and a second application of 5 yields 

5mm = (C,lnX, lnYo)- 

By the symplectic invariance of the F g , we therefore have, order by order in powers of 

9s, 

F g {SMNi) = F 9 (Smm) ■ 

Since our matrix model was engineered to reproduce the Gromov-Witten invariants of Xo, 
we have arrived at 

GW g (X )=F g (C,\nX,\nY ), 

with X and Y obeying the algebraic equation 

H (X,Y ) = 



which coincides with the equation ( 2.3[ ) describing the mirror curve of Xo. 

Given our minimality assumptions on the spectral curve, we have thus derived the BKMP 

conjecture for the fiducial geometry 3£ - 



5.5 The small q limit and the thickening prescription 

The above derivation of the spectral curve for the matrix model is not fully rigorous, as 
we have relied on making minimal assumptions along the way. Although the spectral 
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curve we have found here satisfies all the constraints of section 4.2.2, to prove that it is 
the spectral curve of our matrix model requires a uniqueness result which we currently 
do not have. 

In this section, we provide a heuristic argument that the qualitative behavior of the 
spectral curve and the mirror curve coincide at small q. 

At small q, only very small partitions contribute to the matrix integral. Almost all 



eigenvalues of M, are frozen to the values q a ^ i+d l . By the arguments in section 5.2.1 the 
resolvent Wi(x) hence behaves at small q as 



d 

Wi(x) ~ ^2 ^2 "■ -m-^ + sma11 cut near Q aj ' i+d 

3=0 1=1 



Pictorially, the size of the cuts is shrinking in this limit, replacing the spectral curve by 



its skeleton, see figure 10 



On the other hand, the mirror curve is a priori a tree level quantity, hence does not depend 
on q = e~ 9s . However, recall that we have defined the Kahler parameters Q associated to 
a curve C as 

Q = q J c J . 

The large q limit hence corresponds to the large curve class limit, i.e. the distance between 
the vertices of the pairs of pants out of which the mirror curve is constructed is taken to 
infinity. Just as the spectral curve, the mirror curve thus collapses to its skeleton in the 
q — y limit. 



6 The general BKMP conjecture 

So far, we have obtained the BKMP conjecture only for the fiducial geometry X Q . Studying 
the behavior of the partition function under flop transitions will allow us to extend our 
argument to arbitrary toric geometries. 

6.1 Flop invariance of toric Gromov-Witten invariants 

Under the proper identification of curve classes, Gromov-Witten invariants (at least on 
toric manifolds) are invariant under flops. Assume the toric Calabi-Yau manifolds X and 
X + are related via a flop transition, <f> : X — > X + . In a neighborhood of the flopped 



(—1, —1) curve, the respective toric diagrams are depicted in figure 11 

The 1-cones of S^, corresponding to the toric invariant divisors of X, are not affected by 
the flop, hence can be canonically identified with those of X + . The 2-cones Tj in these 
diagrams correspond to toric invariant 2-cycles Ci, Cf in the geometry. The curve classes 
of X push forward to classes in X + via 

0*([O)]) = -[C O + ], MP& = [<?] + [<%]. (6.1) 
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Figure 10: In the small q limit, only very small partitions contribute to the matrix integral, therefore 
the density of eigenvalues of Mi tends to the flat density (a Dirac comb of equidistant delta functions), 
the non-flat part, which reflects the cuts of the spectral cut, shrinks to zero. 
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Figure 11: X and X + in the vicinity of the (-1,-1) curve. 

— * 

All other curve classes C of X are mapped to their canonical counterparts in j£ + . Under 
appropriate analytic continuation and up to a phase factor (hence the oc in the following 
formula), the following identity then holds 



i.e. 



— * : — * 

Zgw{%~> Qoi Qi, • • • ) Qi, Q) oc Z GW (X , l/Qo, QoQi, . . . , QoQa,, Q) , 



GW g (X, Qo,Qi,-.., Qa, Q) = GW g (X + , 1/Q , Q Q 1 , . . . , Q Q 4 , Q) 



6.2 Proof of flop invariance via mirror symmetry 



Flop invariance of Gromov-Witten invariants upon the identification (6.1) is immediate 



upon invoking mirror symmetry, as (6.1) maps the mirror curve of X to that of X + . The 
proof is a simple computation. 
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Figure 12: X and X + in the vicinity of the (-1,-1) curve. 



Let us introduce the notation to, ti,t$, t^ for the Kahler volume of the curve classes Cj, Cf 
corresponding to the respective 2-cones. In terms of these, we obtain for the mirror curve 
oiX 

Xo + Xl+X2 + ^l e T + d e -n =0j 

Xq Xq 

while the mirror curve of X + is given by 

X±X 2 _T+ x l x 3 -T+ 



Upon invoking £ 3 



x + x ± + x 2 + 



x\X2 —T. 



xo 



Xq 



-e ° + 



0. 



■T2 



e o , we easily verify that the identification (6.1) maps these 



curves and their associated meromorphic 1-forms A into each other. 



6.3 The BKMP conjecture 

Any toric Calabi-Yau manifold X with Kahler moduli Q can be obtained from a sufficiently 
large fiducial geometry (Xq, Qq) upon performing a series of flop transitions and taking 



unwanted Kahler moduli of Xq to oo, see figure 13 for an example. 




Figure 13: Example: We obtain local P 2 from the fiducial geometry with 2x2 boxes by performing 
five flops and then sending the Kahler parameters of the unwanted edges to oo. 

The Kahler moduli of X are related to those of Xo by some relation Q = f(Qo). We have 
just argued that the mirror curves of X Q and X are equal upon this identification, 
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as are the respective Gromov-Witten invariants, 

GW g (X,Q) = GW 9 (X ,Q ). 
Given the BKMP conjecture for the fiducial geometry, 

GW g (X o ,Q o )=F g (S Xo>0o ), 
its validity thus follows for any toric Calabi-Yau manifold: 



GW g (X,Q) = F g (S~a) 



7 Conclusion 

Taking our matrix model from [lj as a starting point and imposing certain minimality 
conditions on the spectral curve, we have thus derived the BKMP conjecture, for closed 
topological strings, for all toric Calabi-Yau manifolds in the large radius limit. As we 
have emphasized throughout, elevating our procedure to a formal proof of the conjecture 
requires a more rigorous derivation of the spectral curve of our matrix model. 

It should also be possible to extend our argument to open Gromov-Witten invariants by 
invoking loop operators, which relate closed to open invariants. In [3], such an operator 
was defined in the matrix model context. An analogous operator should also exist in the 
theory of Gromov-Witten invariants jH] . Establishing the equivalence of these two loop 
operators would allow us to conclude that the W„ s of the spectral curve S% are the open 
Gromov-Witten invariants of X. 

Finally, our treatment of the BKMP conjecture took place at large radius. One should 
study the behavior of the matrix model as one moves away from large radius e.g. to 
orbifold points, and see whether the phase transitions of the topological string are captured 
accurately by the matrix model. Of course, the main tool on the topological string side 
employed in this work, the topological vertex, is no longer applicable in these regions of 
moduli space. 
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A The matrix model 

In this appendix, which is mainly a reprint of section 4 of [1] , we present the matrix model 
which reproduces the topological string partition function on the fiducial geometry 3£o, 
and whose spectral curve we derive in the text. 

SConsider the fiducial geometry 3£ of size (n + 1) x (m + 1), with Kahler parameters 
U,j = a i,j — a i,j+i, r i,j = a i,j+i ~ a i+i,j, an d Si j, as depicted in figure [TJ We write 

Assume that the external representations are fixed to a m+ i = (oio,m+i, a i,m+i, • • • , a n,m+i) 
on the upper line, and a$ = (ao,o, ai,o, • • • , at n ,o) on the lower line. For the most part, we 
will choose these to be trivial. 

We now define the following matrix integral 2 MM ( MM for Matrix Model), 

m+l „ m+l 



hit 1 /• iu-r± p 

Z M M(Q,g s ,dm+i,c%) = A(X(a m+1 )) A(X(a )) ] [ / dM t TT / dR 

m rn 

TT e ^ tr [VxiM-Vti-iWi)] Y[e^ tr t^-i ( M -i)-^W-i)] 

m+l m 

TT e itr(^-M-i)ft "Q 



i=l i=l 

m 



trln/o(A/o) tr ln/ m+ i(M m+ i 



-i) TTgtrln/iCAii) _ / A _-|_n 

i=l 

All matrices are taken of size 

N = (n + l)d. 



d denotes a cut-off on the size of the matrices, on which, as discussed in section 5.2.1 
the partition function depends only non-perturbatively. We have introduced the nota- 
tion 

X{d m+1 ) = diag(X(a m+1 ) i )i=i,...,Af , X{a m+1 ) jd+k = g M°y, m +i) ; 

X(a ) = diag(X(a )i)i=i,...,jv , X(a ) jd+k = q h ^ a ^\ 
for k = 1, . . . , d, j = 0, . . . , n, where 

h i (j)='n-i + d + a. (A.2) 
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Apf) = Yli<j(Xi — Xj) is the Vandermonde determinant. The potentials Vg(x) are given 
by 

n 

V s (X) = -g s J2^(9(q aj /X)) (A.3) 

3=0 

in terms of the g-product 



oo 1 

For z = 1, . . . , m, we have defined 



n=l 



ij-JSL^ It. rf , ~„ 1 -°1.i1 (InQsg " J '')) 2 






(| g(xq l a i,i)g(q a ^/x) 



The denominator of these functions induces simple poles at x — q a ^ i+l for j — 0, . . . , n and 
I G Z. The numerator is chosen such that they satisfy the relation fi(qx) = fi(x). This 
enforces a simple /-dependence of the residues taken at x — q a ^ i+ \ given by a prefactor 
q l - a fact which will be important in the following. These residues are in fact given by 



Res i f l (x)=q a -* +l f jti = -q a ^ +l ] f - 



1 + a.- j — at. j- ,,2 
(In (g J' 1 ro >* )) 



(1)2 e (i + I ^)(l+a J , l -a M )ln 9 ^ ^ 



"; , i " " -*- a o(q a i' i afc .») (1 — q a k,i a j,i\q(q a k,i a j,i\ 



q"l,* 



(A.4) 



where /^ is independent of the integer Z. 
The parameters <% are defined by 



Si — S(j,i-1 + to,i— 1 — s i,«-l — ^_^tk y i + / ik,i-\ ■ (A-5) 



The final equality holds for arbitrary j pQ . 
For % — and i = m + 1, we define 

/ofr) = ^ n — d 



frn+l{x) 



YlUUUix - q h ^) 
1 

n- =0 nti(^-^ (a - +i) ) 



Notice that if the representations c?o or a m+ i are trivial, i.e. hi((Xj t o) = d — i + aj ;0 or 
/ii(<x,, m+ i) = d - i + %, m+1 , we have 

fM-fT ^^' "^ f M-fT ^ xgl ~ aj ' m+1 " d) 
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respectively. The functions f and f m +i have simple poles at x = q hi w>°> (resp. x 
gMoy.m+i)) for I = 1, . . . , d, with residue 

Im = g5So) /o(x) = n jVi nti(? fc,(aw) -9 hi(aj '' o) ) n ¥J (? h,(aw) -? hi( ^ o) ) ' 

11 

fj,m+l;l = Res / m +l(x) 



fyOj.m + l) 



n ^ nf-i(?' ii(aj:m+i) - g fti(a J , .™+ i) ) rL^(<? Maj ' m+i) - g^^.^+i)^ 



The / dependence here is more intricate than above, but this will not play any role since 
the partitions a^o and atj )m +i are kept fixed, and not summed upon. 

The integration domains for the matrices Ri are Hn(R.+ ), i.e. the set of hermitian ma- 
trices having only positive eigenvalues. For the matrices M iy i = 1, . . . , m, the integration 
domains are H^Ti), where 

n 
3=0 

jj t i is defined as a contour which encloses all points of the form g a i> i+N ? and does not 
intersect any contours 7^, (j,i) 7^ (k,l). For this to be possible, we must require that 
the differences a^ — afji be non-integer. The normalized logarithms of two such contours 



are depicted in figure 14 




Figure 14: Two contours surrounding points a + N and b + N, such that a - b £ Z. 



We have defined 

H N (T l ) = {M = UAU\ UeU(N), A = diag(A 1 , . . . , X N ) e T t } , 

i.e. Hfir(Ti) is the set of normal matrices with eigenvalues on Tj. By definition, the 
measure on H^iTi) is (see [24J) 



dM = — A(A) 2 dUdA 



(A.6) 



where dU is the Haar measure on U(N), and dA is the product of the measures for each 
eigenvalue along its integration path. 

The integration domains for the matrices M , M m+ i are Hn(T ), ifjv(r m+ i) respectively, 
where 

n n 

To = (£ i Jfi ) N , r m+1 = (J2 ^m + i) N ■ (A.7) 

3=0 3=0 
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